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The massive star binary fraction in young open clusters I. 
NGC 6231 revisited 
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ABSTRACT 

We present the results of a long-term high-resolution spectroscopy campaign on the O- 
type stars in NGC 6231. We revise the spectral classification and multiplicity of these 
objects and we constrain the fundamental properties of the 0-star population. Almost 
three quarters of the 0-type stars in the cluster are members of a binary system. The 
minimum binary fraction is 0.63, with half the 0-type binaries having an orbital 
period of the order of a few days. The eccentricities of all the short-period binaries 
are revised downward, and henceforth match a normal period-eccentricity distribution. 
The mass-ratio distribution shows a large preference for 0-l-OB binaries, ruling out the 
possibility that, in NGC 6231, the companion of an 0-type star is randomly drawn from 
a standard IMF. Obtained from a complete and homogeneous population of 0-type 
stars, our conclusions provide interesting observational constraints to be confronted 
with the formation and early-evolution theories of O stars. 

Key words: binaries: close - binaries: spectroscopic - binaries: general - stars: 
early- type - open clusters and associations: individual: NGC 6231 ~ open clusters and 
associations: Sco OBI 



1 INTRODUCTION 

Although massive 0-type stars are often so bright that they 
can be studied with small telescopes up to distances of a 
few kilo-parsecs, our understanding of these objects is still 
fragmentary. Even their physical properties remain often ill- 
constrained. The latter however provide key observational 
guidelines to solve one of the most critical astrophysical 
questions at present : how do massive stars form ? In this 
context, we have undertaken a long-term monitoring of the 
O-type star population in a number of young open clusters 
in order to accurately derive and/or confirm their physical 
and orbital properties. This should help to clarify whether, 



one proposed in earlier works, both in terms of the properties 
of the individual objects and those of the 0-star population 
as a whole. 

Located in the core of the Sco OBI association, at about 
1.64 kpc {DM = 11.07±0.04,'Sa na et al.|2006b| |, NGC 6231 
is one of the rich nearby clusters in terms of number of 
hosted O-type stars ([Garcia fc Mermilliod||2001[ |GM01[ ). 
Sana et al. ( 2006b 2007b I identified several hundreds of 



pre-main sequence (PMS) stars, from which they derived 
a cluster age of 2 to 4 Myr, compatible with the evolution- 
ary status of the massive stars in the cluster, and in good 
agreement with previous studies ( Sung et al.||1998 Baume 



as suggested by some authors (e.g. [Penny et al.| 19931 |Garci'a| N aL[1999^ . Using the X-ray properties of these PMS stars to 



I&: Mermilliod|[2"001[ ), a correlation exists between the prop- 
erties of a cluster and those of its member massive star pop- 
ulation. The series of papers initiated here will focus on one 
of the fundamental properties of these objects, their multi- 
plicity, although other properties of the population will be 
discussed whenever allowed by the data set. The present pa- 
per deals with a supposedly well known cluster, NGC 6231. 
Yet, the obtained picture is significantly different from the 
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disentangle the cluster low-mass members from the numer- 
ous field stars along the line of sight, they showed that the 
geometric center of the cluster is located about 30" East of 
the massive binary HD 152248 (Fig.[l]). Using a King profile 
fit, they further constrained the cluster core radius at about 
dc = 3.1' which, at the distance of the cluster, corresponds 
to 1.5 pc. 

Within 15' (= 7.1 pc) around the cluster center (about 
5 times the core radius), one encounters over 90 B-type 
stars, 15 O-type stars and a Wolf-Rayet (WR) system (for 
a complete census, see Sana et al. 2006c I. The spectro- 



scopic binary (SB) fraction of these objects has already 
been investigated by several authors. |Levato fc Morrell| 
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Table 1. Number of radial velocity (RV) measurements obtained for each object with respect to the considered spectral lines. '— ' 
means that the corresponding line has not been measured. 
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1 19831 |LM83| ) derived a minimum fraction of binaries among 
main-sequence early-B/late-0 stars around /min ~ 0.41, 
from which the authors estimated a binary fraction around 
/ ~ 0.57. |Raboud| |T996| obtained a minimum fraction 
/min ~ 0.52 for the B-type stars. More recently, [GMOl] pro- 
posed an even larger fraction of 0.79 for the 0-type stars, 
making NGC 6231 the second cluster in their sample, after 
IC 1805 (/ ~ 0.80), to have an 0-star binary fraction close 
to 80%. In line with the quoted studies, we adopt here the SB 
fraction (/) to be equal to the ratio of the number of binary 
or multiple systems to the total number of objects in the 
considered population. Beyond the SB fraction /, the frac- 
tion of companions per massive star, as defined by|Preibisch| 



et al. (2001 1, will also be briefly discussed in Sect. 



4.5 



In the present paper, based on more than 150 high 
signal-to-noise ratio (SNR) high-resolution spectra, we re- 
visited the current knowledge of all the O-type stars in 
NGC 6231. Our work is organised as follows. The observ- 
ing campaign is briefly described in Sect.p] The individual 
objects are discussed in Sect.[3]while Sect.jljsummarizes the 
properties of the O-type star population seen as a whole. It 
also discusses the implications in terms of constraints on the 
star formation process and on the cluster structure. Finally, 
Sect. [5] briefly summarizes our results. 




Figure 1. O/WR-typc objects in NGC 6231. The circles are 
centered on the geometric center of the cluster and have a re- 
spective radius of 1 X , 2 X , 3 X and 4 X the cluster core radius 
= 3.1' = 1.5 pc. 



2 OBSERVATIONS AND DATA HANDLING 

Our team started to collect data on the O-type stars in 
NGC 6231 about 10 years ago, focusing in the first few years 
on the brightest objects before later extending our survey 
to fainter targets in the cluster. Except for a few spec- 
tra obtained with the Coude Echelle Spectrograph (CES) 
at the ESO Coude Auxiliary Telescope (CAT, La Silla) in 
1998, and with the Bench-Mounted Echelle spectrograph 
(BME) attached to the CTIO 1.5m Ritchey-Chretien Tele- 
scope at Cerro Tololo, most of the present data set was 
acquired between May 1999 and May 2004 with FEROS 
(Fiber-fed Extended Range Optical Spectrograph) succes- 
sively mounted at the ES0-1.5m and ESO/MPG-2.2m tele- 
scopes at La Silla. The data reduction techniques applied 
are identical to those used in previously published analyses 
of the bright short-period binaries in the cluster. We re- 
fer to our previous works on these objects for a complete 
description of the instrumental setups and of the reduc- 



tion techniques (see e.g. [Sana et aL]|2001[ |2003[ |2007a |. 

For each spectrum, Doppler shifts and equivalent widths 
(EWs) of a series of lines (Table [l]) were measured by fit- 
ting Gaussian profiles to the spectral lines. For this purpose, 
we ad opted the effective rest wavelengths fromJConti et al.| 



( 1977 1 below 4800 A and from 



Underhill 



(19941 above that 



limit. The spectral properties of the O-type stars were re- 
derived using the quantitative criteria of |Conti fc Alschuler| 
(19711), IContil (119731), iMathysl (|1988l| arid iMathysl (119891), 



that rely on the EW ratio of diagnostic lines. We adopt the 
usual notations: logW = logiy(A4471) - log VK(A4542), 
logW" = logVy(A4089) - logW^(A4144) and logW" = 
logW^(A4388) + logW(A4686). With respect to the latter 
criterion, based on the product of the EWs (expressed in 
mA), we emphasize that the star has to be single, or that 
the brightness ratio of the two components has to be known 
for this criterion to be applicable. Measured values are re- 
ported in Table [2] Whenever possible, B-type companions 
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Table 2. Diagnostic line ratios and adopted spectral classifications for the studied O-type 
objects. Inappropriate criteria for specific cases are marked as 'n/a'. 



Object 


logVK' 


log W" 


log W" 


Sp. Type 


HD 152076 


0.52 


0.08 - 0.09 


5.3 


09.5 HI 


HD 152200 


0.689 ± 0.048 


0.162 ±0.030 


5.438 ±0.022 


09.7 V 


HD 152233 prim 


-0.337 ± 0.024 


n/a 


n/a 


05.5 ni(f) 


HD 152233 sec 


0.060 ± 0.076 


n/a 


n/a 


07.5 HI/V? 


HD 152234 prim 


0.815 ±0.028 


0.341 ± 0.054 


n/a 


09.7 I 


HD 152234 sec 


~ 0.17 


^ 0.0 


n/a 


08 V 


HD 152247 prim 


~ 0.3 


~ 0.28 


n/a 


09 HI 


HD 152247 sec 


> 0.7 


< 0.0 


n/a 


09.7 V 


HD 152249 


0.39 ±0.02 


0.554 ±0.004 


4.73 ± 0.06 


09 lb ((f)) 


HD 152314 prim 


0.265 ±0.011 


~ 0.27 


n/a 


08.5 V 


HD 152314 sec 


n/a 


n/a 


n/a 


Bl-3? 


HD 326329 


0.49 ± 0.04 


0.03 ±0.03 


5.43 ± 0.02 


09.5 V 


HD 326331 


0.14 ± 0.04 


0.53 ±0.10 


n/a 


08 ni((f)) 


CPD-41°7721 


~ 0.41 


-0.14 - -0.23 


~ 5.5 


09 V 



were also classified by comparing the measured EWs with 



the typical EWs, quoted by Didelon ( 1982 1, for a given spec- 
tral type. 

Table 3 summarizes the derived radial velocities (RVs) 
and provides, object by object, the journal of the observa- 
tions of all the O-type stars, with the exception of the 5 
short-period binaries, that have been the subject of ded- 
icated papers, and of the WR system. References for the 
latter objects are given in Table |4] Tables [l] and 3 show 
that, for most objects, we are able to investigate variability 
on time scales of several hours to several years. 

To further study the possible long-term variability, we 
also compared our measurements with results from earlier 
works. In the past, RV measurements on cluster stars have 



mainly been performed by 


Struve (19441, Hill et al. ( 


1974 


HCB74I, 


Levato & Morrell ( 


1983 LM83I, Perry etal. ( 


1990 


PHYB90 


1 and Garcia & Mermilliod (2001 GMOll. 


Note 



however that these authors report (mean) RVs based on dif- 
ferent lines, so that direct comparison of the quoted RVs 
from one author to the other, or with the current work, 
should be considered as indicative only. The most recently 



published photometric studies are from Perry et al. ( 1991 



PHC91I, Balona & Laney (1995 



(1997 RCB97I, Sung et al. (1998 



BL95I, Raboud et al. 



SBL98I and Baume 



et al.| ( |1999[ |BVF99[ ). Sung (2005, private communication) 
also obtained UBV(RI)c photometry over a 40' x 40' field 
of view down to ~ 23 (for a brief description of the data, 
Sana et al.|2006b[ ). 



Investigating the multiplicity of massive stars is a dif- 
ficult task, not only because of the numerous observational 
biases, but also because of the large parameter space that 
must be searched. In the present approach, we adopt either 
of the following two criteria as a definite proof that an ob- 
ject is a binary system: (i) the presence of Keplerian (and 
periodic) RV variations, or (ii) the detection of a composite 
(SB2) and variable spectrum. 

In the first case, we are able to compute an orbital solu- 
tion, thus constraining most of the orbital parameters. In the 
second case, the detection of the secondary spectrum also 
provides valuable information. Indeed, the spectral types of 
both components offer a reasonable estimate of the mass ra- 



tio and the detected RV variations yield a first, albeit rough, 
idea of the time scale of the orbital period. 

As described below, some objects in our sample show 
RV variations that we have not been able to link with a 
Keplerian motion, nor have we been able to detect the sec- 
ondary signature. While some of these objects are poten- 
tial binary candidates, their multiplicity nevertheless awaits 
more definite confirmation. Combined with the fact that we 
are actually observing the entire massive star population in 
the cluster and not a sub-sample of it, the present approach 
offers the advantage of providing a firm lower limit on the bi- 
nary fraction in NGC 6231, which is unaffected by statistical 
errors. 



3 THE O-TYPE STAR POPULATION IN 
NGC 6231 

As mentioned earlier, 15 O-type objects are found within 
5 core radii from the cluster center (Fig. [T]). In the present 
section, we revisit the properties of each of these objects 
individually. 



3.1 Short period binaries 

Detailed analyses of the 5 short period binaries were pre- 
sented in a series of papers since 2001. Table [4] summarizes 
the main physical and orbital properties of these objects. We 
note that they significantly differ from the results quoted in 
the contemporaneous work of |GM0l| Implications for the 
general properties of the NGC 6231 O-star population, as 
well as constraints on our understanding of the formation 
and dynamical evolution of these objects, will be discussed 
in Sect.H 



3.2 Long period binaries 

3.2.1 HD 152234 

HD 152234 is reported in the literature as a BOIab radial ve- 
locity variable star ( Levato fc Malaroda|1980 K Levato et al 
(19881 proposed a very preliminary period of 27.25 days, 
while 



GMOl I claimed to have found the companion signature 
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Table 4. Selected orbital and physical parameters of the short period binaries in NGC 6231. 



Object 



P (d) 



mi/m2 



Sp. Type 



Reference 



CPD-41°7742 
HD 152219 
HD 152248'' 
HD 152218 
CPD-41°7733 
WR79(HD 152270) 



2.44070 

4.24032 

4.81602 

5.60391 

5.681504 

8.8911 



0.027" ±0.006 1.803 ±0.015 09 V + B1.5 V 



0.082 ±0.011 
0.133 ±0.005 
0.259 ± 0.006 
0.0 (fixed) 
0.0 (fixed) 



2.530 ±0.023 
1.012 ±0.011 
1.319 ±0.014 
2.640 ±0.012 
2.7 



09 III + Bl-2 V/III 
07 III (f) + 07.5 III (f) 
09 IV + 09.7 V 
08.5 V + B3 
WC7 + 06V 



Sana et al.| | j2003 
Bana et al. (2006a| 
Bana ct al. ( 201 
Sana ct al. ( 200i 
Bana ct al. ( 2007a| 
E uehrs| | |1997| >;f 




° We note that |Bouzid et al.lpOOS^ reported a circular orbit for CPD -41°7742. 

° The orbit has been recomputed using the new version of the Liege orbital solution package |Sana et aT]|2006a| for a 
better handling of the errors. 
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4460 4480 5860 5880 

X (A) X (A) 



4536 4544 

X (A) 



Figure 2. HD 152234: HelAA4471, 5876 (left) and Henl4542 
(right) line profiles at three different epochs. The orbital shift of 
the primary is clearly seen. The secondary component is clearly 
visible, blueshifted, in the top spectrum and redshifted in the 
bottom one. Vertical dashed lines indicate the adopted rest wave- 
lengths. 
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Figure 3. HD 152234: Hell4471 RVs folded against two values 
of the orbital period. Filled (resp. open) symbols indicate primary 
(resp. secondary) RVs. 



3.2.2 HD 152233 



HD 152233 is known as an 06111(f) star (Walborn 1972 



Levato & Malaroda 19801 that displays small RV shifts 
(e.g. Struve 19441. Based on data from the hterature 
along with four additional observations, |GM01| proposed 
a short period of 4.15 d. According to their Table 2, they 
measured the same RV in two spectra separated by one day. 
Consequently, they adopted an extremely large eccentricity 
e = 0.57, thus making HD 152233 a potential cornerstone 
for the study of the dynamical evolution of O-type binaries. 



and derived a period twice as long (P — 54.6 d). Their com- 
puted orbital solution is slightly eccentric (c = 0.183±0.093) 
and they estimated Ki — 46.5 ± 2.8 kms~^ and 
K2 = 211.5 ± 6.1 kms"^(r.m.s. = 13.5 kms~^) . 

Between 1999 and 2004, we have obtained 23 FEROS 
spectra and 5 BME spectra. From our observations, HD 
152234 most probably consists of an 09.7 supergiant pri- 
mary and an 08 main sequence secondary. The secondary 
signature is faint (Fig.[2| and is only seen at very few 
phases. This partly results from the lower SNR of some 
of the spectra in our data set. Our current data definitely 
rule out the period of 54.6 d proposed by |GM01| and 
point towards a preliminary value close to 126 d (Fig. [3|. 
Although the orbital properties await further constraints, 
we nevertheless confirm that HD 152234 is a long period 
O-type binary. 



In May 2004, we observed the system twice a night 
for six consecutive nights, thus well sampling the proposed 
4.15 d period. However, we definitely ruled the latter out eis 
we only detected a slow trend in the system RVs (Fig. |4|. 
The current data set points to a period of a few hundred 
days but awaits further constraints. Our May 2003 data 
revealed for the first time the signature of the secondary 
component in the blue wing of the He I and He 11 lines (Fig. 
[5|. Despite the fact that the lines are only marginaly disen- 
tangled, we estimated individual spectral types of 05.5 and 
07. 5 for the two components. The Henl4686 line displays 
a variable P-cygni profile. Although the current data set 
does not allow us to separate the relative contributions 
of the two components to the Henl4686 profile, its EW 
is compatible with both stars being giants. Finally, we 
also detected Nlll lines that are usually not present in a 
typical 05.5 HI spectrum. These lines are moving in the 
opposite direction compared to the primary lines. Clearly 
this is also the signature of the companion, a less evolved 
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Figure 4. HD 152233: He I AA 4471 and 5876 line profiles as ob- 
served during a 6-night run in May 2004. HJD at mid-exposures 
are quoted at right-hand in the format HJD— 2 450 000. The ver- 
tical dashed lines represent the rest wavelength. No significant 
difli'erence in the position of the absorption lines is observed from 
one spectrum to the other, ruling out the 4.15 d period derived 
by |GM01| 



O star of spectral type 07.5 or later. While more numerous 
disentangled spectra are needed to better constrain the 
nature of the secondary, HD 152233 is definitely a long 
period O-l-O binary. 



3.2.3 HD 152247 

At 11' N of HD 152248, HD 152247 has clearly suffered from 
a lack of attention in the previous works. Mostly classified 
so far as 09.5 HI (|Houck|| 
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Figure 5. HD 152233: He 1 14471 and He II 14542 line profiles 
as observed during different observing runs from 1999 to 2004. 
Clear shifts are seen and the profile asymmetry in May 2003 
(HJD~2 452 784) suggests that the secondary is a late O-type 
star. The layout of the figure is similar to the one of Fig. |4] 



stant in the V band ( Balona|1983 1 with a range in the pub- 
lished values smaller than 0.05 mag. Owing to the absence 
of the star in the usual photometric works on the cluster, we 
adopted the magnitude of Diplas & Savage (1994): V = 7.18. 
The first RVs were obtained by Struve^ (4944^ who quoted 
values around —16 kms~^(rms~10 kms~ ). No other mea- 
surement was obtained until Raboud ( 1996 1 suggested that 



the star was displaying RV shifts. He indeed reported 5 mea- 
surements spread over two observing epochs (see the webda 
databas^. The first set was acquired over three consecutive 
nights and indicated values around -f20(±3) kms^^. The 
next two measurements were obtained about two years later, 
again on two consecutive nights, and yielded values c l ose to 
-23(±2) kms"^ More recently, [Stickland fc Lloyd| (|2001 1 



Mathys[1988 Penny et al 



p56l [Feinstein fc Ferrer] p68} 



19961, the star seems to be con- 



reported two additional observations from the International 
Ultr aviolet Explore r (lUE) obtained about 10 years be- 
fore [Raboud] ([^96| and indicating RV ' 39(±1) kms"\ 

Though the star clearly showed RV variations, the time scale 
over which these occur remained very poorly constrained. 

In our data, clear RV shifts are noticeable from year-to- 
year (Figs. [6] and [7|. The two extreme RV values in our set 
are approximately -40 (1999, 2004) and -1-20 kms"^ (2006). 
In May 2004, we obtained six spectra over six consecutive 
nights but the HD 152247 spectrum remained almost un- 



http: / / www.univie.ac.at / webda 
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Figure 6. HD 152247: He I AA 4026, 4471 and Hell AA 4542, 4686 line profiles as observed during different runs from 1999 to 2006. The 
spectral signature of the companion is clearly seen in the 2006 spectrum (at HJD~2 453 836). The layout of the figure is similar to the 
one of Fig. [4] 
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Figure 7. RV variations of the three long period binaries with ill-constrained periods: HD 152233 (Hell5876), HD 152247 (Hell5876) 
and HD 152314 (Helll4686). 



changed. Similarly no significant shift was observed between 
the 2002 March and April data. In April 2006, we unam- 
biguously detected for the first time the signature of the 
secondary star. Values for log W' and log W" point towards 
the primary being an 09 III star. In addition, the secondary 
displays a clear Helll4686 line in absorption and is thus an 
0-type star as well. Because of the blending with neighbour- 
ing primary lines, the disentangling of the secondary clas- 
sification lines is difficult. To correct for the primary line, 
we subtracted a spectrum where the two components were 
deblended from one where they were blended, taking into ac- 
count the measured Doppler shift. In principle, the primary 
signature should be removed, leaving a positive and negative 
image of the secondary spectrum shifted from one another. 



Using this rough method, we estimate that the secondary is 
probably an 09.7 dwarf. 

Again, more disentangled spectra are needed to confirm 
the classification of the secondary and to constrain the 
orbital period, which is probably of the order of years. Still, 
the present results indicate that HD 152247 is an additional 
O+O binary in NGC 6231. 



3.24 HD 152314 

At 4.'2 E-NE of the cluster cor e, HD 152314 is a late O-t ype 
object. With V = 7.92 - 8.12, [Feinstein fc Ferrer] |T968l) re- 
ported that the star displays photometric variability. |RCB97l 
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Figure 8. HD 152314: HelAA4026, 4471, Hell AA 4542, 4686 and Hell5876 lines acquired during five different observing runs. The 
spectral signature of the companion is clearly seen in the May 2003 spectrum (at HJD~2 452 784). The layout of the figure is similar to 
the one of Fig. |4] 



further found long-term variability on a time-scale of ~ 18 
years. These variations display different amplitudes in the 
different bands of the Geneva system: the star is becom- 
ing bluer when it brightens. |RCB97| also suggested an ad- 
ditional, shorter time-scale variability of the order of 1.5 



months. RV measurements by |Struve (19441, PHYB90 and 
|LM83| do not present any particular variation and the lat- 
ter authors reported the star as RV constant. More recently, 
[GMOT] obtained three additional spectra: two separated by 4 
days, another acquired about two years later. Though with 
some uncertainty on the RV measurements, they reported 
to have detected the signature of a blue-shifted secondary 
component on the first spectrum of their series, thus suggest- 
ing variability on a time-scale of days. The spectral types 



quoted in the literature are OS. 5 III (Levato & Malaroda 
1980|, 09 III (|Schild et al.|1969[[PeFry et al.|1990[ ) or 09 V 



i jMorgan et al.||1953a[ ). The remark by |RCB97| that the ob- 
ject's colour is changing with time could also indicate a slight 
change in its spectral type. 

We have acquired a few FEROS spectra of HD 152314 
since May 2000. In May 2004, we acquired 9 spectra over a 
6-night run. The night-to-night variability is very limited 
but we observed significant changes from one year to 
the other (Figs. [7] and |8|. In May 2003, we detected the 
presence of a secondary companion in most of the He I and 
Balmer lines. The Ha profile presents mixed absorption 



and emission and evolves from an inverse P-Cygni profile 
to a normal P-Cygni through the years, a behaviour that 
could be compatible with an absorption component moving 
on top of a slightly broader emission. The spectral type 
changes from 09 to OS. 5 between 2000 and 2004, but 
this could be due to an effect of line blending. The best 
spectral type estimate is probably obtained in 2003, when 
the two components were deblended, which indicates an 
OS. 5 HI primary. The Hell lines seem to belong solely to 
the primary. The secondary component is thus probably 
a B-type star. However, its signature remains undetected 
in the metallic lines, rendering an accurate spectral clas- 
sification very difficult. Because of the strong secondary 
signature in the He I lines, its classification could be located 
around the He I maximum, i.e. at spectral-types B1-B3. 
This classification remains however uncertain. Finally, with 
V = 7.749 and B - = 0.525, we obtained Mv = -4.95 
for the system, signific antly fainter than a typical OS. 5 HI 
star {M^^-^™ = -5.32, [Martins et al.|2005| . This situation 
is reminiscent of the one observed for other O-type binaries 
in the cluster. We thus assume that the primary is rather 
a dwarf and we adopt, as our best classification, 08.5 V 
-I- Bl-3 V. This object definitively deserves follow-up 
observations covering time scales from weeks to years. 
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Figure 9. HD 152200: Left panel. Evolution of the H/3 and He 1 14921 lines during six consecutive nights in May 2004. The layout of 
the figure is similar to the one of Fig.|4] Right panel. He 1 14471, He II 14686 and He 1 15876 hues obtained in 2000 (HJD = 2 451 669.900, 
sohd Une) and 2004 (HJD = 2 453 135.676, dashed line). 



3.3 Variable stars 



3.3.1 HD 152200 



HD 152200 is a late O-type star located about 3.'6 W-SW 
of HD 152248. It is quoted as variable by |Feinstein fc Fer-| 
Q ( |1968[ ), with V ~ 8.35 - 8.42. Photoelectric measure- 
ments reported in other works since the 1960's are all in the 



range 8.41-8.44, except the measurement from Oja ( 1986 
V 



37). More recent CCD photometry yielded values be- 
tween 8.31 ( |Baume et al.|1999[) and 8.4 2 ( |BL95[ ). Regarding 
the spectral types, Sch ild et al.| ( |1969 ') classifie d HD 152200 
as 9.5 HI , [Houk (1978[) a s Q9 HI, whereas L ovato fc Malar- 



oda| ( |1980[ ) and |PHYB90| respec tively preferred 09.5 V(n) 
and BO IV classifications. |LM83| reported three RV measure- 
ments on three consecutive nights, which indicate a veloc- 
ity change of about 70 kms~^ f rom one n ight to the other. 
The Ca-K corrected velocities of |PHYB90l show a 30 kms"^ 
range on three measurements spread over 10 days. Finally 
|GM01| claimed to have detected the secondary component 
signature on 2 out of their 5 spectra. The system should ap- 
parently have a period of a few days as the lines are reported 
as blended one night and as separated by ~ 180 kms~^ the 
following night. 

From May 1999 until May 2004, we collected 16 FEROS 
spectra spread over different time scales. The optical spec- 
trum of HD 152200 is characterized by the Balmer, He I and 
Hen lines. The few metallic lines typical of late O-type 
objects are also clearly seen. The line profiles are variable 
(Fig. [9| but no obvious signature of a companion could be 
found. The observed variations display a peak-to-peak am- 
plitude of 50kms~^ and seem to have a period close to 5 
days. The H/3 line behaviour (Fig.[9]l is clearly different from 
what was reported by .GMOl,. The maximum night-to-night 
shift is about 30 kms~ (to be compared with the 70 kms^^ 
shift observed by LM83I. From our observations, the cen- 



tral value around which the variations occur seems to have 
shifted by about 10 to 15 km s"'^ between the years 2002 and 



2004. It could however be an apparent effect resulting from 
a too sparse sampling of the possible ~ 5-day period. We 



used the Lafler & Kinman ( 1965 \ and the Heck et al. ( 1985 1 



period search techniques and two aliases were found around 
P ~ 4.44 and 5.56 days. However, both period values yield a 
large scatter of the RV measurements plotted against a hy- 
pothetical phase and it seems difficult to interpret the RV 
variations in terms of a Keplerian motion. Still, the ampli- 
tude of the RV variations compared to the ones of the other 



O-stars in the cluster (Fig. 10 \ pinpoints HD 152200 as a 
good binary candidate. 

The measured EWs correspond to an 09.7 star with 
spectral type 09.5 at 1-a. Conti's criterion {W") clearly in- 
dicates a giant class while Mathys's criterion {W'") rather 
points towards class V. We computed the absolute visual 
magnitude to be about My = —4.08 which is typical of a 
main sequence 09.7 star. Though binarity would provide the 
most straightforward explanation to the observed RV vari- 
ations, we have not been able to satisfactorily fit the data, 
nor have we been able to detect the companion's signature. 
More observations are needed to confirm or invalidate the 
multiplicity status of this object. 



3.3.2 HD 152249 

HD 152249 is one of the brightest members of the clus- 
ter. Its visual magnitude, as quoted since the 1960's, is 
in the range 6.43-6.51, with the notable exception of |Fe-| 



instein & Ferrer ( 1968 1 who reported the star as variable 
Quoted spectral types are 09- 
who listed 07- 



with V = 6.34 - 6.50 
BO la/b/ab, with the exception of 



RCB97 



09 V ! HD 152249 is adopted by [Walborn fc Fitzpatrick] 
( |1990[ ) as an OC9.5 lab standard. Several authors indicated 
RV variations. Neubauer (19301 obtained a ARV range of 
46 kms~^ on 3 plates and classified it as a spectroscopic 
binary. |Struve| ( |1944[ ) and [Perry et al.| ( |1990[ ) respectively 
obtained 15 and 13 measurements on a time-span of 16 and 



10 H. Sana et al. 



75 days. Both data sets display a similar range as the one 



obtained by Neubauer ( 1930 1. However these do not present 



the smooth variations that could be expected for a binary 
with a period larger than a couple of days. Later on, |Gar-| 
many et al. ( 1980 1 and Levato et al. ( 1988 1 obtained ad- 



ditional measurements that do not indicate night-to-night 
variability but, separated by about 5 years, they differ by 
20 kms"^ 

Between 1999 and 2004, we collected 34 spectra. The 
peak-to-peak dispersion over the 6 years is 20 kms~^ 
and the night-to-night variability has an amplitude of 
~ 10 kms~^. We consider the star to present slight RV 
changes but clearly, we cannot associate these with a binary 
nature. The measured EWs unambiguously lead to an 09 I 
spectral type, thus slightly different from the [Walborn] 
fc Fitzpatrlck] ( |1990| ) classification. With V = 6.44 and 
B — V = 0.19, we obtained Mv = —6.1, quite typical of lb 
supergiants rather than lab ( |Lang||T992| | . With the presence 
of the N III AA 4634-4641 lines in weak emission, we finally 
adopt an 09 lb ((f)) classification. 
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3.3.3 HD 326331 = CPD-imU 



With V ~ 7.5 (|SBL98|, HD 326331 is an O-type star 
located at 3' East of the cluster center. Together with 
CPD-4r7744B {V ~ 10), it forms a visual pair separated 
by 7.3". If physically related, their revolution period should 



be over 200 000 years (Mason et al. 19981. Quoted spectral 



types range from 07 III (Levato & Malaroda 1980 Lev- 
|ato fc MorreII|p83| ) to Q9 III ( |HCB74[ ), but only give the 
composite classification. The system was reported as photo- 
metrically variable by |PHC91[ and |HCB74| mentioned an 
amplitude of 0.3 mag. If we ignore the value V — 7.71 



from Feinstein & Ferrer (19681, then all the published pho- 



tometric data since the 1960's indicate a maximum am- 
plitude of 0.15, with an average value of 7.52 ± 0.08. RV 



measurements were obtained by Struve (19441, PHYB90 



|HCB74| and |LM83| The |HCB74] data consisted of 13 spec 
tra obtained over 7 consecutive nights. Based only on the 
11 14069 and Si iv 14089 lines, they reported smooth veloc- 
ity variations between -1-94 and —124 kms~^ (mean error 
of 13 kms~^) and then up again. They thus suggested that 
the star is a binary with a period P > 7 d. They how- 
ever cautioned that "the spectrum is very difficult to mea- 
sure as the lines are weak". Finally, using one additional 
measurement obtained about 7 years later, [LM83| published 
a first SBl orbital solution with P = 6.24208 ± 0.00015, 
e = 0.12±0.09 and Ki = 75±8 kms"^ The obtained resid- 
uals (~17 kms~^) are however very large compared to other 
orbital solutions published in the same paper. 

We have sparsely observed HD 326331 since 1998. We 
obtained 4 spectra of the He 1 14471 region with the CES 
in May 1998, which revealed very broad lines (> lOA). We 
then successively acquired FEROS spectra in May 1999 (5), 
in April 2002 (2) and in May 2004 (6). These were usu- 
ally separated by one night to cover the proposed 6.4-day 
orbital period. The spectrum of HD 326331 presents the 
clear signature of the usual Balmer, He I and Hen lines. 
Both He II 14686 and Ha profiles show mixed absorption and 
emission. A faint emission component could also be present 
in the red wing of He 1 15876. The lines present a fiat bot- 



Figure 10. Mean RV (airv) and l-<j dispersion ((Trv) for the 10 O- 
type stars/binaries described in Sect. |3.2"] and [33| Circles indicate 
the definite binaries in our sample. This graph is based on the 
He 115876 line measurements only to preserve the homogeneity of 
the data set and to avoid biases induced by combining different 
lines for different objects, ctrv is arbitrarily set to for the 2 
objects with only 1 RV measurement. 



tom and their profiles are therefore clearly not Gaussian 



(Fig. 11 1. The broadness of the lines definitely suggests that 
HD 326331 is a rapid rotator and wo used the Fourier trans- 
form of the line profile ( (Gray 2005, and references therein) 
to estimate a rotation velocity wsin(i) — 310kms~^. Then, 
assuming that the shape of the lines is fully determined 
by the rotational broadening, we used a rotation profile 
with vsm{i) = 310 km s"^ to fit the HelAA4471, 5876 and 
Helll4542 lines. The obtained RVs are reported in Table 3 
and show 1-a range of about 3 to 9kms~^. 

Though the peak-to-peak dispersion seen in our mea- 
surements is about 20 kms~^, we never obtained such highly 
positive RVs as the ones reported by HCB74 Still, Fig. pT] 
shows that the spectrum of the star displays obvious profile 
changes. Though their nature is unclear, it seems difficult to 
relate them to the presence of a secondary companion and 
we consider them as intrinsic to the star. 

To estimate the spectral type, we measured the 
EWs of the usual classification lines. For this peculiar 
object, we integrated the different line profiles rather than 
adopting the EWs of the fitted profile. Measured values 
point towards an 08 I star. However, with the Hen 14686 
and Ha lines showing mixed absorption and emission, 
with N III AA 4634-41 weakly in emission and with clearly 
delineated He 1 14388, the observed spectrum is not that of 
a supergiant. Furthermore, rapid rotation is rather unlikely 
for a supergiant. We thus prefer the III((f)) classification. 
With V = 7.50 and B - V = 0.18, the visual absolute 
magnitude of the object is about My = —5.2, which again 
rules out a supergiant classification. According to |Howarth| 
& Prinjal (|1989| and iHumphreys & McElroyl d 19841), this 
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Figure 11. HD 326331: He I AA 4471, 5876 line profiles as ob- 
served from 1999 to 2004 on various consecutive nights. The lay- 
out of the figure is similar to the one of Fig. |4] 



value is rather typical for giant stars. While determining 
the nature of the observed LPVs requires more follow-up, 
we can however rule out the 6 or 7 day period previously 
claimed by |LM83| and, until proven to be otherwise, we 
consider the star as single. 



3.4 Presumably single stars 

34.1 HD 152076 

Somewhat offset from the cluster core, HD 152076 had so 
far received less attention than the other bright stars. Some 
confusion seems to exist about its magnitude and indeed, 
cross-identifiers reported in the SIMBAD database are clearly 
erroneous. For example, SIMBAD associates HD 152076 with 



It is obvious from a comparison with neighbouring objects in 
DSS images that the V = 10.85 magnitude reported in SIM- 
BAD is wrong. We adopt in the following V — 8.471 (Sung 
2005, private communication) which is close to previous de- 
t erminations of ^ = 8.5 ( IPHYB90|), 8.48 (|PHC91[ ), 8.47 
(Schild ct al. 1969), 8.50 Feinstein & Ferrer 1968) 



8.46 



Bok et al.,,1966J or 8.48 ( |Heske fc Wendker„1984J Quoted 



spectral types are in the range BO V ( Morgan et al.|1953a|b I 
- BO/1 HI ( [Schild et al.|pCT| |Houk| |1978[). A couple of RV 



measurements were also performed by Struve ( 1944 1, Wilson 
(^953]) and[PHYB90]who obtained a Ca- K corrected RV of 



Balona 



( 1983 1 reported this ob- 



about —25 to —30 kms" 
ject to be constant in the Johnson B filter on a time scale 
of 5 hours. 

We acquired one snapshot FEROS spectrum of 
HD 152076 that reveals relatively narrow lines. The Hell 
signature is faint, but clearly seen, and spectral crite- 
ria indicate an 09.5 V/III star. With V = 8.471 and 
B — V = 0.240, we obtained My = —4.4, an intermediate 
value between typical absolute magnitudes of 09.5 dwarfs 
and giants. HD 152076 is quite offset from the cluster core 
and thus might not belong to NGC 6231. Instead the star 
could belong to the Sco OB 1 association, located at the 
same distance and with a similar age as NGC 6231. Though 
its luminosity does not agree perfectly with the spectral 
type deduced from the spectroscopic criteria, HD 152076 is 
unlikely to be a background object. Because of the narrow 
lines, suggesting a lower gravity, and the relatively strong 
metallic spectrum, we opted for an 09.5 HI classification. 
RV measurements performed on the classification lines 
plus Mg II 14481 give RV = -30.7 ± 2.7 kms~\ in good 
agreement with previous measurements. We thus consider 
HD 152076 to most probably be a single star, though it is 
obviously difficult to discard the possibility of a very long 
period binary. 



3.4.2 HD 326329 = CPD~il°7735 

HD 326329 is located in the core of the cluster. With 
V ~ 8.8, it is relatively bright and has thus been observed 
for a long time. Derived spectral classifications oscillate be- 
tween 09 and BO and agree on the main sequence luminos- 
ity class. A relatively large scatter in the reported visual 
magnitudes led some authors to consider this object as vari- 
able. Focusing on CCD data published since the 1990's, we 
find V in the range 8.71 ( |RCB97| ) 



1.81 ( |BL95| ). [LM8 3 and 

|PHYB90| respectively obtained 3 and 2 spectra of HD 326329 
and found consiste nt RVs (respectively RV — — 30(±2) and 
— 30(±5) kms~^). GMOl obtained three additional spec- 
tra spread over 6 days. They derived RVs from —34 to 
—60 kms~^, thus largely scattered and quite different from 
the cluster systemic velocity, but did not comment on these 
facts. 

We obtained two FEROS spectra in May 2000 and 
six additional spectra in May 2004. Within our data set, 
HD 326329 does not present any significant RV change. For 
example, we obtained RV = -20.3 ± 1.5 kms"^ for the 
Hell4388 line, RV = -18.0 ± 1.9 kms"^ for the Hell5876 
line and RV = -16.5 ± 2.2 kms"^ for the He 1 17065 line. 
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Figure 12. Left panel: Cumulative distribution of the number of O-type objects (solid line and open symbols) and O-type binaries 
(dashed line and filled symbols) as a function of the distance d from the cluster geometric center. Middle panel: Evolution of the binary 
fraction with d. Right panel: Integrated O-star number density (cr = N/nd'^) versus d. Symbols and lines have the same meanings as 
in the left panel. In all three panels, vertical dotted lines indicate multiples of the core radius dc = 3.1' = 1.5 pc. 



The spectrum of HD 326329 displays the clear signature 
of an 09.5 V star, with the 09 type at 1-a. Finally, with 
V = 8.76 and B - y = 0.17, we obtained My = -3.8, which 
is slightly too faint for an 09.5 V star at the cluster distance. 

Though no significant RV shift is found from our data 
alone, we note that the average RVs reported by different 
authors are only marginally compatible with our data. This 
could suggest a long period modulation but only a long 
term monitoring of this object could bring an answer to 
this question. 



strengths of the Si ill 14552 and Si iv 14089 lines are very 
similar. Comparing with the spectral atlas of |Walborn fc| 
Fitzpatrick (19901, we finally adopt a Bl V class. With 
V = 8.709 (resp. V; = 9.892) and {B - V) = 0.168 (resp. 
{B - V), = 0.178), the visual magnitude of CPD -41°7721 
(resp. CPD-41°7721s) is about A/y = -5.2 (resp. 
My.b = —2.6). Both stars are a few tenths of a magnitude 
too faint if located in the cluster core. Alternatively, slightly 
later spectral types would yield a much better agreement. 



34.3 CPD~4r7721 

Together with CPD-41°7721s {V = 8.90), CPD-41°7721 
{V = 8.72) is a visual double star with components 
separated by ~5'.'8 and corresponding to SBL 350 and 
351 or B VF 12 and 27. Two spectra were obtained by 
PHYB90I indicating RV = -23 and -28 km s"^ (mean error 



-10 kms" 



GMOl 



also acquired two spectra of the bright 
component and quoted RV — —26.8 and —33.0 kms^^. 
They assigned it a spectral type 09.5 V. 

In May 2004, we obtained one FEROS spectrum of 
CPD-41°7721 and another of CPD -41°7721s. Each of 
the two spectra show a single spectral signature. Conse- 
quently, we consider the two stars as apparently single. 
CPD-41°7721 is an 09 V star with an average RV 
of —25.5 kms~^, although we observed a large scatter 
(from —36 to —13 kms~^) according to the various lines 
measured. With no He 11 14542 line and a faint He 11 14686 
line, CPD— 41°7721s is clearly an early B-type star. The 



4 DISCUSSION 

4.1 Minimum SB fraction 

As mentioned earlier, and taking into account WR 79, 16 ob- 
jects containing at least one O-type star are located within 
15' from the cluster center. As shown in Figs.[l]and 12 these 
stars are unevenly spread across the field of view. Most of 
the objects are located in the vicinity of the cluster center: 
ten of them lie within one cluster core radius (dc ~ 3.1'), 
while 13 of them are to be found within 1.5 x dc. Among 
these 16 objects, 6 are close SB2 binaries with periods be- 
low 10 days, 4 are longer period systems, 3 stars display line 
profile variations (LPVs) and 3 stars show constant RV and 
constant line profile. 

Given the evidence presented in Sect. |3j no doubt re- 
mains about the multiple nature of the 6 short-period and 
of the 4 longer period systems, though the orbital proper- 
ties of the latter ones await tighter constraints. This yields 
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a firm lower limit on the fraction of binaries in NGC 6231 of 
/min ~ 0.63. 

In Sect. [3j we further noted that other objects were 
presenting hints of binarity, among which HD 152200 and 
HD 326329 are probably the two best candidates. If the 
multiplicity of these objects is confirmed, this would signif- 
icantly increase the binary fraction compared to the lower 
limit mentioned above. Yet we consider that, so far, there 
is not enough observational evidence to support such an as- 
sertion. 



4.2 Statistical significance of the results 

While the previous paragraph established that in NGC 6231 
the true O-type binary fraction / is located between 0.63 
and 1.00, the present results have implications beyond the 
sole case of NGC 6231 and allow us to put some interesting 
constraints on the parameters of the underlying distribution 
of single vs. binary stars. 

Let us adopt a binomial distribution B(n, p), with n 
the number of trials and p the probability of realisation, to 
describe the current experiment (i.e. counting multiple stars 
in a sample of n objects). Given the current realisation of 
the experiment (10 multiple stars among 16 objects), one 
can compute that the probability of realisation p of the par- 
ent distribution has to be larger (at the 0.01 significance 
level) than 0.37. In the hypothesis that some binaries have 
remained undetected despite our efforts, the parent distri- 
bution parameter p would be drawn towards larger values, 
so that the lower limit p > 0.37 still holds at the 1% signifi- 
cance level (or better). 



4.3 Orbital parameters 

The orbital parameters that we have derived for the short 
period binaries in NGC 6231 (see Table [4| also reveal a sig- 
nificantly different picture than previously accepted. In par- 
ticular, while GMOl proposed eccentricities up to 0.6 for the 
short period binaries, we discard these values as resulting ei- 
ther from erroneous Julian dates that pull the eccentricities 
towards larger values (see e.g. Sana et al.|2003 1 or from ill- 
constrained periods biasing the orbital fit. From our results, 
none of the short period binary stars has an eccentricity 
larger than 0.3. As a consequence, the NGC 6231 binaries 
are now in line with the period-eccentricity distribution of 



other O-type binary objects (Fig. 13 1. This has significant 
implications on the dynamical evolution of the O-type stars 
as no other binaries with a high eccentricity and a very short 
period are known so far. 

We also note that the short-period systems are signif- 
icantly more abundant than longer period systems, even 
compared to Opik's law which states that the distribution 
of the periods (or equivalently, given Kepler's third law, of 
the separations) should be flat in the logarithmic space. In- 
deed, this would translate itself into a linear dependence of 
the cumulative distribution of the systems, expressed ver- 
sus the logarithm of the period while Fig. |14| shows that, in 
NGC 6231, this is not the case. Even if all the non-binaries 
in the cluster were actually undetected long-period binaries, 
this fact would remain, so that observational biases cannot 
be invoked. It is however difficult to decide whether this is a 
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Figure 13. Period - eccentricity distribution for the O-type SB 
systems listed in the 9th Spectroscopic Binary catalogue ( |Pour-| 
|baix et al.||2004l l. Filled symbols indicate the location of the 
NGC 6231 systems. 
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Figure 14. Cumulative number of binary systems, in NGC 6231, 
as a function of the natural logarithm of the orbital period. Period 
values above ln(P) > 4 are very preliminary, though the order of 
magnitude should be correct. 



signature of the formation process or of the early evolution 
history. 
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4.4 Mass-ratios and companion IMF 

Not only is the present study based on a complete and ho- 
mogeneous population of 0-type stars, but we emphasize 
that the signature of all the secondary companions could be 
unveiled in the detected binaries, yielding an unprecedented 
view of the binary population parameters. As a matter of 
fact, all the 10 detected O-type binaries have an OB com- 
panion, suggesting luminosity ratios [L\/ L2) between 1 and 
10, and mass ratios (m\/ra2) between 1 and 3. 

Of course the present study is biased towards the detec- 
tion of relatively massive companions. Nonetheless, taking 
advantage of the completeness of the population, one can 
already make the following strong statement. Even if the 6 
remaining presumably single stars are undetected binaries 
(low mass or faint companions, pole-on view, etc.), a min- 
imum of two thirds of the companions have a mass above 
5Mq. Even considering the possibility of hierarchical mul- 
tiple systems, the present data set definitely rules out the 
possibility that the companions in NGC 6231 are randomly 
drawn from a Salpeter initial- mass function (IMF). Indeed, 
if this were the case, one would expect, for each O+O sys- 
tem, to encounter 9, 13, 19, 22, over 100 and a few hundreds 
of O type binaries with a companion of spectral type B, A, F, 
G, K and M respectively ( Weidner fc Kroupa|2008" I . Consid- 
ering that, in the present work, one would probably detect O 
and B companions only, and assuming a 100% binary frac- 
tion, a rule of thumb would predict a binary detection rate 
of 2% only, thus about 0.3 binary given the 16-star sample. 
Despite the limited size of the sample, the null hypothesis 
that the companions of an 0-star are randomly drawn from 
a Salpeter IMF can be rejected at the 5-cr levej^ 



4.5 Visual companions 



HD 152249 have a physical faint companion with putative 



revolution periods between 150 and 4900 yr (Mason et al 



The speckle interferometric campaign of Mason et al. ( 1998 1 



included the brightest objects {V < 8) of the cluster and re- 
vealed that HD 152233, HD 152234, HD 152248 and, possi- 
bly, HD15224^have one or several visual companions with 
a sub-arcsec separation. Beside these results, direct imag- 
ing indicates that CPD-41°7721s is located at 5'.'8 from 
CPD -41°7721, that HD 326329 has a visual companion at 
7'.'3 and that HD 152076 is not distant by more than 8" from 
two other stars. Only HD 152200 remains seemingly alone in 
a 20" radius. So far, however, it is almost impossible to state 
whether these visual pairs are physically linked, or arise by 
line-of-sight coincidence in such a crowded field. 
Using only stars brighter than V = 17, one can estimate a 
density of about 10 stars per arcmin^ at the cluster center 
and about half this value at 1.5 x dc ( [Sana et al.||2006bl ). 
Rough statistics thus suggest that companions separated 
by less than 1 and 2" respectively (according to whether 
they are located in the very core or in the outer part of the 
cluster) are unlikely to arise by chance and are thus likely 
bound (at the 0.01 significance level). According to this sim- 
plistic criterion, only HD 152234, HD 152248 and, possibly, 



p98l). 



If verified, this would indicate that at least one third 
of the SB systems are hierarchical triple systems, yielding a 
number of companions per massive star close to 0.8. Such a 
value remains still far below the 1.4 companion fraction de- 



rived by Preibisch et al. (2001 1 for the O- and B-type stars 



in the Orion Nebula Cluster. As a matter of comparison, 
either a significant fraction of close and unresolved compan- 
ions, with a physical separation roughly in the range 2-200 
A.U., are still missing, or one would need to consider all the 
known visual companions up to 10" to reach a similar rate 
for NGC 6231. 



4.6 Physical constraints for massive star 
formation 

Because of its spectroscopic approach, the present study is 
clearly biased toward the detection of SB systems with an 
OB star companion and only provides a lower limit on the 
true binary fraction. Even considering the work of [Mason] 
et al. (19981, about 2 orders of magnitude in semi-major 
axis (i.e. between ~ 2 and 200 A.U.) are still not sampled. 
Investigations of this range would require 10 m/s RV ac- 
curacy or milli-arcsec spatial resolution, both of which, as 
far as O-type stars are concerned, challenge the current in- 
strumentation limits. Still the present study provides strong 
constraints for massive star formation; 

- At least 60% of the O-type objects in NGC 6231 are SB 
systems. Over 70% of any O-type stars in the cluster are 
found to belong to an 0-1-0 or an 0+B system. 

- If the properties of the O star population in NGC 6231 
are representative of the O star population in general, the 
expected binary fraction of the parent distribution should, 
at the 0.01 significance level, be larger than 0.37. 

- Concerning the 10 SB binaries in our sample, half of 
them have an O-type companion, the other half having a 
BO-3 companion. Therefore, in all the detected SB binaries 
in the present sample, the luminosity ratio is about or less 
than 10 to 1 while the mass ratio is about or less than 3 to 
1. 

- Even in the extreme case of a 100% binary content, it 
is very unlikely that the companion of an O-type star would 
be randomly drawn from a standard IMF. 

- As a corollary, the previous statement shows that the 
massive star formation and/or early evolution tend to pro- 
duce a large number of double objects with components hav- 
ing masses of the same order of magnitude. 

- 60% of the detected SB systems (thus about one third of 
the O-type objects) are tight binaries with periods between 
2.4 and 9 days, and with zero or small eccentricities. Com- 
pared to Opik's law, short period systems are overabundant, 
even in the extreme case of / = 1.0. 

- No significant difference in the spatial distribution of the 



binary fraction can be found across the cluster (Fig. 12 1 



^ We note that using a standard Kroupa IMF instead of a 
Sa lpeter IMF do es not significantly change the present result. 
^ [Mason et al.| suggested that HD 152249 could have a binary 
companion at a separation of 0'.'06, but they emphasized that 
this result needs to be confirmed. 



5 CONCLUSIONS 

Based on a set of high-resolution high SNR spectra of the 
O-type objects in NGC 6231, covering various time scales 
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from a few hours to several years, we have revised their 
physical and orbital status, paying particular attention to 
their possible multiplicity. The present fraction of binaries in 
NGC 6231 (around 63%) is significantly reduced compared 
to the value of 79% obtained by |GM01| 

We note that the revision of the binary fraction in 
NGC 6231 puts the cluster in a less extreme position com- 
pared to other rich (in terms of the number of 0-type stars) 



open clusters. Considering also that De Becker et al. ( 2006 1 



could not reproduce the high binary fraction (/ ~ 0.80) 
found by |GM01| for the cluster IC 1805 either, the anti- 
correlation proposed by |GM01| between the cluster density 
and its massive star binary fraction, remains to be con- 
firmed. 

While detailed studied of other open clusters are 
needed, the present approach already offers the possibility 
to quantitatively test the null hypothesis that the O-type 
SB populations of different clusters are drawn from the same 
distribution. Indeed, the rejection of this hypothesis is a pre- 
requisite to the suggested link between the binary fraction 
and the cluster properties. 

Beyond the binary fraction, we also conclude that, in 
NGC 6231, the massive stars are not randomly paired from 
a standard underlying IMF. Instead, there is a strong bias 
towards the formation of O-I-OB systems. Among the latter, 
about half of them are close binaries with a period of a 
few days. Finally, none of the short period binaries has an 
eccentricity larger than 0.3. Compared to the periods and 
eccentricities found by |GM01[ this has strong implications 
on the dynamical evolution of such systems. 

The present results outline the limitations of previ- 
ous estimates of the properties of the O-type population in 
young open clusters. It also emphasizes the need for exten- 
sive studies to accurately constrain, among other properties, 
the binary fraction of various stellar populations and, par- 
ticularly, of early-type stars whose formation and evolution 
are still not firmly understood. 
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